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LINEAR  GEOMETRY 

We  consider  first  the  electrostatic  problem  of  a  rectangular  beam  in  a 
straight  rectangul ar  box  ( cgs>  esu). 
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By  integration  over  x,  from  w  -  h  to  w  +  h  ,  and  by  integration  over 

X  A  X 

y1  from  -hy  to  y  using  V>  and  from  y  to  hy  using  V<>  the  potential 
function  within  the  (uniform)  beam  becomes 
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wherein  it  will  be  recognized  that  sin  y  is  equal  to 
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-1 
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etc. 


for  n  =  1 


etc. 


As  a  check,  we  note  that 

V2V  -  -16  p  V  l  sin  Si  sin 
i—t  n  2 


n'hX  .in 

2w  "ST 
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0  for  0  <  x  <  w  -  h  and  for  w  +  h  <  x  <  2w. 

X  X 


In  particular,  if  hx  =  w, 


V^V  =  -16  p 


_sin£*7sin  3S+Fs1n  5S  +  •••] 


=  -16  p  j  =  — 4irp  (for  0  <  x  <  2w). 


[cf .  B.O.  Pierce,  Eq.  808.] 


If  we  form 
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When  self  magnetic  forces  are  also  taken  into  account  in  this  straight-pipe 
situation,  correspondingly  the  gradient  of  effective  field  becomes  multiplied 

O 

by  the  factor  1  -  e  ,  where  6  =  v/c. 
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TOROIDAL  GEOMETRY 


This  toroidal  case  has  been  analysed  by  Dr.  Lloyd  Smith.  For  the 
electrostatic  problem  we  now  employ,  in  our  present  notation,  radial  functions 


F0(";r)  -  a0  [jQ(knr)  Yfl(kna)  -  YQ(knr)  JQ(kna)] 

*  ‘o  Yo  -  T«fIo,B)t)  J“  f“<n)t)]  • 

where  the  kn,  or  zQ(n)  =  kp  "R,  are  such  as  to  make  FQ(n;r)  vanish  for  r  =  b  and 
oQ  is  an  arbitrarily  selected  factor  (e.£. ,  aQ  =  1). 

The  Green's  function  then  may  be  written: 
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Integration  over  the  location  rp  yp  of  the  source  involves  first  the 
following 
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Followed  by  an  integration  over  r,  (from  R  -  h  to  R  +  h  ). 

1  XX 
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For  this  latter  integration  we  follow  the  suggestion  of  Dr.  Smith  in  that  we 

take  p(r,  )  =  p  — —  ,  whereupon  we  obtain 
1  rl 
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Finally,  then  —  for  this  toroidal  electrostatic  problem  —  if  we  form 
,  we  obtain 
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Correspondingly,  for  the  magnetic  problem,  we  obtain  a  similar  result  —  save 
that 

(i)  a  function  Fj_  is  employed,  in  which  the  Bessel  functions  are  of 
order  1, 

(ii)  the  zeros  are  the  quantities  z^(n)  for  which  F^(n;b)  =  0, 


(iii)  a  factor  is  to  be  appended  for  the  contribution  to  the  force 
gradient: 
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These  terms  [(2a)  and  (2b)]  combined  give  the  expression  proposed  by 

Dr.  Smith,  if  we  identify  our  pR  with  his  ep  R. 

o 

In  the  work  that  follows,  we  shall  identify 
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Following  what  we  believe  to  have  been  the  suggestion  of  Dr.  Smith,  we  have 
undertaken  to  evaluate  the  integrals  (and  the  subsequent  sums)that  appear  in 
the  above  expressions  numerically.* 


*4 — order  integration  by  use  of  the  extended  Simpson's  Rule  (Abramovitz  and 
Stegun,  §25.4.6) 
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We  remark  that  TERM  1  and  TERM  2  will  be  rather  close  to  one  another 
numerically  —  and  indeed,  when  w/R  <<  1,  will  each  be  rather  close*  to  the 
quantity 
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exhibited  by  expression  (1)  [p.  2]  divided  by  4irp. 


♦Thus  if  Fq  (or  F^)  be  approximated  by 
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if  we  identify  x  =  r  -  (R  -  w) 
=  r  -  a 
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— thus  providing  agreement,  in  this  limit,  with  expression  (4).  Note  that 
terms  for  which  n  is  even  will  be  expected  to  be  small  [since  the  factor 
sin  n  ir/2  in  the  above  approximation  vanishes  for  n  even].  Also,  moreover, 
if  hx  =  wand  x  =  w(r  =  R),  this  approximation  suggests  that 


Such  approximate  behavior  has  proven  of  use  in  checking  some  aspects  of 
operation  of  the  computer  program  CPRTS. 


COMMENT: 

We  recall  that  for  a  straight  beam  pipe  the  electric  and  magnetic  self 
forces  are  expected  virtually  to  cancel,  through  the  action  of  the 

I  2 

factor  -y  =1-8  . 

y 

In  the  toroidal  geometry,  on  the  other  hand,  the  electric  and  magnetic 

2 

forces  may  differ  slightly  in  magnitude  even  in  the  limit  e  >  1. 

To  compare  the  relative  sizes  of  such  effects,  it  may  be  helpful, 

therefore,  to  think  of  a  quantity  analogous  to  that  we  form  from  the 

y 

toroidal  solutions  as 

TERM  1  -  TERM  2 
Avg.  of  TERM  1  and  TERM  2 

—  1*  ®.*  > 

— - - -  .  (TERM  1  -  TERM  2). 

TERM  1  +  TERM  2 
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THE  PROGRAM 


The  program  CPRTS  (Library  LASLETT)  first  computes  the  first  100  "zeros" 

[zQ(n)  and  z-|(n)  of  FQ(n,b)  and  F-j(n,b)]  by  successive  application 

of  Newton's  method  —  employing  in  this  connection  the  Sub-routine 

(Mrs.  Barbara  (Harold)  Levine)  LJBJY  for  Bessel  functions  and  (optionally) 
s  t 

their  1 — derivatives.  (The  chamber  radii,  a  and  b,  are  entered  for  this 
purpose.)  Optionally  printed,  in  tabular  form,  are  the  zQ(n)  and  then  the 
z-j(n),  each  accompanied  by  the  value  to  which  FQ  or  F-j  has  thereby 
been  reduced  to  virtually  zero. 

Next,  with  hx  entered,  the  program  then  computes  (and  prints) 


R+h 
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R-h 


Fo<n;rl >drl 


J  F1(n;r1)dr1 
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RHOfll  =  - ^ 


RH01  = 
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I  I 


r  [F0(n;r  )]  ‘dr 


/D  i  o  i 

r'[F1(n;r  )]^dr 


[Integration  steps  are  taken  to  be  inversely  proportional  to  n.] 

Actually  the  numerators  and  denominators  of  each  of  these  quotients  are 
calculated  in  two  ways:  First,  directly;  second  by  subtracting  some  simple 
form  from  the  integrand  and  then  supplementing  the  integral  of  the  result  by 
the  analytic  integral  of  the  simple  form  that  had  been  subtracted.  The 
ratios  (RHO0  and  RH01)  computed  by  the  "direct"  and  by  the  "difference" 
procedure  are  each  printed.  One  then  has  the  option  of  selecting  whether 
one  will  retain  for  future  use  results  (quotients)  obtained  by  the  direct  or 
by  the  difference  procedure  and  such  values  are  then  printed,  once  again,  as 
a  summary.  [For  a  given  number  of  integration  steps  the  results  obtained  by 
the  difference  procedure  are  believed  to  be  somewhat  more  accurate.] 
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By  entering  hc,  h^,  and  the  field-point  radius  (rg),  the  program 
is  then  in  a  position  to  compute  the  successive  terms  (n  =  1,  2,  3,  ...)  and 
to  print  the  successive  cumulative  sums  of  the  expressions  TERM  1  and  TERM  2. 

Upon  the  completion  of  such  a  tabulation,  the  program  can  be  directed  to 
accept  entries  of  new  values  of  hc,  h^,  and  rg  —  or  alternatively  to 
return  to  earlier  stages  of  the  program. 

RESULTS 

Results  for  several  different  parameters  (including  rD)  are  summarized 

D 

in  Tables  that  follow. 
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a  =  990.00 


b  =  1010.00 

hc  =  5.00 


w  =  10.00 
h  =  9.90 

A 

hy  =  4.50 


rB  1 

TERM  1 

TERM  2 

TERM  1  -  TERM  2 
AVERAGE 

1000.00 

_ ! 

0.89152  8556 

0.89151  9000 

1.0719  x  10-5 
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Data  Tabulated  p.ll 
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a  =  87.50 


b  =  112.50 


w  =  12.50 


hc  =  10.00 


h  =  8.40 

A 

hy  =  2.00 


rB 

TERM  1 

TERM  2 

TERM  1  -  TERM  2 

AVERAGE 

92.80 

0.65204  9938 

0.65143  3302 

'  9.46136  x  10“4 

95.20 

0.82374  9845 

0.82296  2369 

9.56422  x  10~4 

97.60 

0.86558  6391 

0.86472  6339 

9.94100  x  10"4 

100.00 

0.86185  8622 

i 

1  0.86100  3476 

f 

9.92704  x  10"4 

102.40 

0.82865  0396 

|  0.82786  4456 

9.48908  x  10~4 

104.80 

0.75552  2275 

I  0.75486  2845 

f. 

8.73183  x  10-4 

107.20  I 

0.57508  3712 

|  0.57460  6194 

8.30690  x  10"4 
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a  =  187.50 


b  =  212.50 


w  =  12.50 


hx  =  12.46 

h  =  10.00  h  =  2.00 

c  y 


rB 

TERM  1 

TERM  2 

TERM  1  -  TERM  2 
AVERAGE 

189.32 

0.49962  2554 

0.49954  6695 

1.51844  x  10“4 

192.88 

0.81219  1130 

0.81201  3801 

2.18358  x  10"4 

196.44 

0.88358  9404 

0.88336  2741 

2.56558  x  10"4 

200.00 

0.88990  1674 

0.88966  5564 

2.65357  x  10-4 

203.56 

0.85477  4512 

0.85456  3070 

2.47396  x  10"4 

207.12 

0.76014  5318 

0.75999  0786 

2.03313  x  10"4 

210.68 

j  0.45297  1393 

* 

0.45290  9272 

1.37151  x  10-4 

a  =  187.50 

rB 

b  =  212.50 

hc  .  10.00 

TERM  1 

w  =  i; 

hx=  * 

"y- 2 

[ 

TERM  2 

2.50  (as  above) 

1.40 

!.00  (as  above) 

TERM  1  -  TERM  2 
AVERAGE 

192.80 

0.63065  1758 

0.63050  80  7  2 

2.27863  x  10-4 

195.20 

0.80521  8479 

0.80503  1029 

2.32821  x  10"4 

197.60 

0.85568  1330 

0.85547  1844 

2.44848  x  10“4 

200.00 

0.86149  2848 

0.86127  9728 

2.47415  x  10"4 

202.40 

0.83723  5836 

..... 

0.83703  5575 

2.39222  x  10~4 

204.80 

0.77118  1650 

0.77101  0104 

2.22470  x  10"4 

207.20 

0.59231  7299 

0.59219  0831 

i 

1.13537  x  10~4 
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a  =  387.50  b  =  412.50  w  =  12.50 

h  =  8.40 

A 

hc  =  10.00  hy  =  2.00 


rB 

t 

TERM  1 

TERM  2 

.  -  -  -  . 

TERM  1  -  TERM  2 
AVERAGE 

392.80 

0.62054  1001 

' 

0.62050  6263 

5.5982  x  10~5 

395.20 

0.79634  6686 

0.796  30  0  904 

5.7492  x  10"5 

397.60 

0.85090  4066 

0.85085  2326 

2.0808  x  10-5 

400.00 

0.86140  1722 

f  '  " 

0.86134  8483 

i 

6.1807  x  10"  5 

402.40 

0.84168  4122 

i 

1  0.84163  3534 

[ . . . 

6.0104  x  10~5 

404.80 

0.77933  7854  i 

:  ! 

|  0.77929  4057 

5.6199  x  10“5 

407.20 

j 

0.60139  2252  | 

:  j 

|  0.60135  9661 

5.4194  x  10' 5 
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